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Abstract 



We present arguments for the existence of five-dimensional rotating black holes with equal magnitude 
angular momenta in Einstein-Gauss-Bonnet theory with negative cosmological constant. These solutions 
posses a regular horizon of spherical topology and approach asymptotically an Anti-de Sitter spacetime 
background. We discuss the general properties of these solutions and, using an adapted counterterm 
prescription, we compute their entropy and conserved charges. 

1 Introduction 

In five dimensions, the most general theory of gravity leading to second order field equations for the metric 
is the so called Einstein-Gauss-Bonnet (EGB) theory, which contains quadratic powers of the curvature. 
The Gauss-Bonnet term appears as the first curvature stringy correction to general relativity [U [5] , when 
assuming that the tension of a string is large as compared to the energy scale of other variables. 

The study of black holes with higher derivative curvature in Anti-de Sitter (AdS) spaces has been con- 
sidered by many authors in the recent years. Static AdS black hole solutions in EGB gravity are known in 
closed form, presenting a number of interesting features (see e.g. [3], [1], [S], [Hj and the references therein). 
It is of interest to generalize these solutions by including the effects of rotation. This problem has been 
considered recently in [7] within a perturbative approach. The authors of [7] discussed some properties of 
a particular set of asymptotically AdS^ (d > 4) rotating solutions in EGB theory with one nonvanishing 
angular momentum (where the rotation parameter appears as a small quantity), the effects of an U(l) field 
being also included. 

The main purpose of this work is to present numerical evidence for the existence of a different class of 
rotating solutions in d = 4 + 1 EGB theory with negative cosmological constant, approaching asymptotically 
an AdS spacetime background. These solutions are found within a nonperturbative approach, by directly 
solving the EGB equations with suitable boundary conditions. They posses a regular horizon of spherical 
topology and have two equal magnitude angular momenta. This leads to a system of coupled nonlinear 
ordinary differential equations (ODEs), which are solved numerically. The same approach has been employed 
recently to construct Einstein-Maxwell rotating black hole solutions in higher dimensions [51 [S] . 

2 The general formalism 

2.1 The action and boundary counterterms 

We consider the EGB action with a negative cosmological constant A = 




(2.1) 



1 



where R is the Ricci scalar and 

LgB= -4.R^.R'"' +R^uarR'"'^\ (2.2) 

is the Gauss-Bonnet term. The variation of the action (|2.ip with respect to the metric tensor results in the 
equations of the model 

1 a 

R^i^ - ^R9t.u + ^9t.u + -^H^^ = , (2.3) 

where ^ 

Hfj,i, = 2{Rfji,aKTRj"^^ — "^Rp.puaR'^'^ — '^R^aR'^ u + RR^v) — -^LoBgfj.i' ■ (2.4) 



For a well-defined variational principle, one has to supplement the action (|2.1|) with the Gibbons-Hawking 
surface term 

1 



dM 



and its counterpart for the Gauss-Bonnet gravity [2] 



^""""^ = - T7^ / d'^^ (J - 2Gab^°') , (2.6) 

where jab is the induced metric on the boundary, K is the trace of the extrinsic curvature of the boundary. 
Gab is the Einstein tensor of the metric jah and J is the trace of the tensor 

Jab = ^{2KKacK'b + K^dK^'^Kab - 2KacK'''Kdb - K^Kab) ■ (2.7) 

To compute the conserved charges of the asymptotically AdS solutions in EGB gravity, we use the 
approach proposed by Balasubramanian and Kraus in [10]. This technique was inspired by AdS/CFT 
correspondence and consists in adding suitable countcrterms let to the action of the theory in order to ensure 
the finiteness of the boundary stress tensor Tab ~ ■^/= derived by the quasilocal energy definition [TI] . 

Therefore we supplement the general action (which contains the surface terms for Einstein and Gauss- 
Bonnet gravity) with the following boundary counterterm 



O^^tj^ JdM ^ 



(2.8) 



where R is the curvature scalar associated with the induced metric 7 (see also [12] for previous work on 
boundary conterterm technique in EGB gravity, applied to non-rotating solutions). The consistency of the 
procedure fixes the expressior0 of Ci, C2: 



ci = ' ^ \ C2 = ^ ^ 3 ^ ■ (2.9) 



Varying the total action with respect to the boundary metric 7ofc, we find the divergence- free boundary 
stress-tensor 

Tab = -^^{Kab~ Kjab + Ci'yab + C2Gab+'^{Qab~^Qlab)^ ' (^-lO) 

where [13] 

Qab - 2KKacK^ - 2KacK'''Kdb + KabiK.dK^'' - K^) + 2KRab + RKab - 2K'''R,adb - ^RacK^ , (2.11) 



^ As a — > 0, one recovers the known expression in Einstein gravity, ci -> -3/e + a/4e^ + 0{a)^ , C2 ^ e/2 + 3a/8i + 0(a)^ 
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with Rabcch Rafc denoting the Riemann and Ricci tensors of the boundary metric. 

Provided the boundary geometry has an isometry generated by a KiUing vector , a conserved charge 

f d^^S' e'Ty (2.12) 

can be associated with a closed surface S PJ]. Physically, this means that a collection of observers on 
the hypersurface whose metric is 7 all observe the same value of 0^ provided this surface has an isometry 
generated by ^. 

2.2 The metric ansatz and known limits 

While the general EGB rotating black holes would possess two independent angular momenta and a more 
general topology of the event horizon, we restrict here to configurations with two equal magnitude angular 
momenta and a spherical horizon topology. The suitable metric ansat^l reads [8| 

ds^ = — - + g(r)d9^ + h(r) sin^ 6 (d(p - w{r)dtf + hlr) cos^ 9 Id^ - w(r)dt? (2.13) 
/('') 

+ (.g(r) - h{r)) sin^ 6'cos2 e{dLp - d^jf - h{r)dt^ 

where G [0, 7r/2], {f.^i) G [0, 27r], r and t being the radial and time coordinates. For such solutions, 
the isometry group is enhanced from R x U{\Y to R x U{2), where R denotes the time translation. This 
symmetry enhancement allows us to deal with ODEs (in what follows, we fix the metric gauge by taking 
9{r)=r'^)- 

For the metric ansatz ()2.13p . the EGB field equations (|2.3p present two well known exact solutions. The 
first one corresponds to the generalization jSj of the static Schwarzschild-AdS solution with a Gauss-Bonnet 
ternll 



f{r)=b{r) = l+'-\^l-^l + 2a(j^-^^y g{r) ^ h{r) ^ r\ w{r) ^ . (2.14) 

The AdSs generalization [16l [l^ of the Myers-Perry rotating black holes [18] with equal magnitude angular 
momenta is found for a = (no Gauss-Bonnet term) and has 

/W = 1 + ^ - ^ + Hr) =r'l^l + — ^ j , wir) = -5^, g{r) = r^ b{r) = ^,(2-15) 

where AI and a are two constants related to the solution's mass-energy and angular momentum, while 

s = 1 - aye. 

3 Black hole properties 

We are interested in black hole solutions with an horizon located at a constant value of the radial coordinate 
r = r/i > 0. Restricting to nonextremal solutions, the following expansion holds near the event horizon: 

fir) = /i(r - r,,) + 0(r - r^)2, hir) = K + 0{r - rh),b{r) = bi{r - m) + 0{r - r^f Mr) = Wh + 0{r - r^)(3.1) 

The event horizon parameters r^j, /i, 61, and (with (/i, 61, h^) > 0) are related in a complicated 
way to the global charges of the solutions. 



^ EGB rotating topological black hole with zero scalar curvature of the event horizon are known in closed form (see e.g. |14j . 
I15| and references there). However, they are found for a different metric ansatz and present rather different properties. 

^Note that the EGB gravity presents two kind of static black hole solutions, which are classified into the plus and the minus 

branches, (r) = fe(-i-) (r) = 1 -|- ^1 ± ^^l + 2a — jj^^ ■ In this paper we shall restrict to the minus branch solutions, 

which present a well defined Einstein gravity limit. 
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The metric functions have the foUowing asymptotic behaviour in terms of the constantsQ /2, 62 and 



W4: 



/ = l + -(l-\/l-— ) + lt + 0(l/O, b=l + r^ ^- + -± + Oil/r% (3.2) 




1 _ /1 2a 



+ 0(l/r6), u;(r) = ^ + 0(l/r«) 



One can see that, similar to the static case, the parameter a must obey a < £ /2, beyond which the theory 
is undefined. For these asymptotics, the effective cosmological constant is Ae// = A(l + y^l — ^')/2. 

The Killing vector x — d/dt + ^^djdip + Vl^d/dil) is orthogonal to and null on the horizon. For the 
solutions within the ansatz (|2.13p . the event horizon's angular velocities are all equal, = fJi^ = Uh- The 
Hawking temperature as found by computing the surface gravity is 



Th = (3.3) 
47r 

Another quantity of interest is the area Ah of the rotating black hole horizon 

Ah - y^r^^a, (3.4) 

where V3 = 27r^ denotes the area of the unit three dimensional sphere. 

These rotating solutions present also an ergoregion inside of which the observers cannot remain stationary, 
and will move in the direction of rotation. The ergoregion is the region bounded by the event horizon, located 
at r = r/i and the stationary limit surface, or the ergosurface. The Killing vector d/dt becomes null on the 
ergosurface , i.e. gu — —b{r) + h{r)w{r)'^ = 0. The ergosurface does not interesect the horizon. 

3.1 The global charges and entropy of solutions 

The global charges of these solutions are computed by using the counterterm formalisnlfl presented in Section 
2. The computation of the boundary stress-tensor Tab is straightforward and we find the following expression 
for the components of interest here 



rjit _ ^ 

" SttG' 



hO(l/r''), (3.5) 



2a 



The mass-energy E of solutions is the charge associated with the Killing vector d/dt, 

^ = £;(0) ^^(c)^ where (3.6) 



SttG 2 V ' ^ 16 1^ P \l 

where E^'^'^ represents the Casimir term |10) in EGB gravity, presenting a nontrivial a dependence^ (this term 
appears also in the static limit (|2.14p '). These black holes have also two equal magnitude angular momenta 
Jip = J-4i = J, with 



^The MPAdSs solution ifzTsj l has /2 = 2M(b? /l'^ - 1), = -2M, Wi = 2Ma. For the static Schwarzschild-AdS-Gauss- 
Bonnet solution 112.141 1 one finds /2 = ^2 = — m/^ 1 — U14 = 0. 

different computation of the mass and angular momentum of EGB solutions was also reported namely in 1191 . 1201 , 1211 . 

°In the small a limit, one finds E^") = 37r£2/32G - 21a7r/64G -I- 0(a)2. Note that the first order correction to the mass of 
pure global AdSs does not depend on the value of the cosmological constant. 
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representing the charges associated with the KiUing vectors d/d(p, d/dtp as computed according to (|2.12p . 
Also, in what follows it is important to use the observation that one can write 



1 

2p 



J- A. 

bh dr 



ll^hw'l -r- 



+ a(/-4+^) 



(3.8) 



^GB) 




r^hww' - 5') + a f (/ - 4 + 4 + 4^)^' + (4 - ^ - f)hww' + rfhw" 



The gravitational thermodynamics of the EGB black holes can be formulated via the path integral approach 
[22l I23j. However, while the static vacuum Lorentzian solutions (|2.14[) extremize also the Euclidean action 
as the analytic continuation in time has no effect at the level of the equations of motion, this is not the 
case of the rotating configurations discussed in this paper. In this case it is not possible to find directly real 
solutions on the Euclidean section by Wick rotating t ^ it the Lorentzian configuration^ In view of this 
difficulty one has to resort to an alternative, quasi-Euclidean approach as described in [24]. The idea is to 
regard the action / used in the computation of the partition function as a functional over complex metrics 
that are obtained from the real, stationary, Lorentzian metrics by using a transformation that mimics the 
effect of the Wick rotation t it. In this approach, the values of the extensive variables of the complex 
metric that extremize the path integral are the same as the values of these variables corresponding to the 
initial Lorentzian metric^. 

When computing the classical bulk action evaluated on the equations of motion, one replaces the R — 
2A + jLgb volume term with 2{Rl + jH^) and make use of (|3.8p to express it as a difference of two boundary 
integrals. A straightforward calculation using the asymptotic expansion (|3.2p shows that the divergencies of 
the boundary integral at infinity, together with the contributions from and , are regularized by 

Jet- As a result, by using also the first set of relations in (|3.8p . one finds the finite expression of the classical 
action 



1 



4G^V7^ 



(/2-4Ma/1-^ + ^(1 



6 a 



P 



2a , 



(3.9) 



4 / 2a 



Upon application of the Gibbs-Duhem relation to the partition function, one finds the entropy S — (3{E 
2ujhJ) — IcU which is the sum of one quarter of the event horizon area plus a Gauss-Bonnet correction 



^3 2 / — 

S = Sa + Sgb, with 5*0 = — r^V^fti 



Sgb = a—y/hh{A- -j). 



(3.10) 



In the static limit, the known expression S 
solutions is 5 = V^Ah /^G. 



^(r^ + 3ar/j) is recovered, while the entropy of the MPAdSs 



4 Numerical results 

The EGB equations (|2.3p lead to a system of four coupled second order ODEs for the metric functionf0 
/(r), 6(r), h{r) and 'w{r). We are interested in solutions of these equations presenting the asymptotic 
expansion p.ip . (|3.2p . 

^ Even if one could accompany the Wick rotation with various other analytical continuations of the parameters describing 
the solution {e.g. a ^ ia for a MPAdSs black hole), given the numerical nature of the configurations in this paper, there is no 
assurance that the modified metric functions will also be solutions of the field equations in Euclidean signature. Instead, one 
has to solve directly the EGB field equations for a metric ansatz with Euclidean signature. 

®Note also that not all closed form solutions with Lorentzian signature present reasonable Euclidean counterparts, in which 
case one is forced again to consider a 'quasi-Euclidean' approach. The d = 5 asymptotically flat rotating black ring solutions 
provides an interesting example in this sense 1251 . 

^These equations are extremely complicated (each of them containing around flfty terms) and we shall not present them 
here. 
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Figure 1. A typical rotating solution with Wh = 1-8 (solid line) is plotted together with a static solution [wh = 0, 
dashed line) for ru = 0.5, a = 0.5, = 20. 
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Figure 2. Two rotating solutions with rn = 0.5, Wh = 0.5, = 20 are plotted for a = 0.1 (dashed line) and 
a = 1 (solid line). 

In order to construct numerical solutions, the constants {a, A) have to be fixed. Then the sohition is 
further specified by the event horizon and the angular velocity at the horizon w{rh) (or cquivalcntly, the 
angular momentum J through the parameter W4). 

The complete classification of the solutions in the space of parameters is a considerable task that is not 
aimed in this paper. Instead, by taking the arbitrary value (."^ = 20, we analyzed in detail a few particular 
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Figure 3. The parameters fi,bi,hh at the event horizon are plotted together with the parameters /2,62,W4 in 



the asymptotic expansion at infinity, as a function of a for solutions with rh ~ 0.5, 



20, Wh = 0.5 (Figure 3a). In 



Figure 3b we plot the Hawking temperature, the mass-energy iS'"', the Casimir term E^'^\ the angular momentum 
J and the entropies So,Sgb for these solutions. Here and in Figure 4b iS'"', iS''^' and J are plotted in units with 
Vs/SttG — 1, while we set V3/4G = 1 in the expression of S. 



classes of solutions, which hopefully would reflect all relevant properties of the general pattern. However, we 
have found nontrivial rotating black hole solutions for other values of the cosmological constant, in particular 
for A = and for A > 0. 

Also, since the Gauss-Bonnet term in (j2.ip has to be considered as a correction to the Einstein-Hilbert 
action, we report here the results for positive values of a in the interval < a < 1 (however, solutions with 
£^ = 20 and larger values of a exist as well) . 
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Figure 4. The parameters fi,bi,hh at the event horizon are plotted together with the parameters /2, 62,^4 in 
the asymptotic expansion at infinity, as a function of the event horizon radius rn for solutions with a = 0.5, £^ = 
20, Wh ~ 0.5 (Figure 4a). In Figure 4b we plot the Hawking temperature, the mass-energy i?'"', the angular momen- 
tum J, the entropy of the solution in Einstein gravity 5*0 and the Gauss-Bonnet correction Sgb for these solutions. 

In the absence of a closed form solution, we relied on a numerical methods to solve the equations. The 
numerical methods here are similar to those used in literature to find other numerical black hole solutions 
vi^ith equal magnitude angular momenta [8l[9]. We take units such that G = 1, and employ a collocation 
method for boundary-value ordinary differential equations, equipped with an adaptive mesh selection pro- 
cedure [55]. Typical mesh sizes include 10'^ — 10** points. The solutions have a typical relative accuracy of 
10-8. 

In constructing rotating EGB-AdS black holes, we make use of the existence of the closed form solutions 
p.l4p and (|2.15p . and employ them as starting configurations, increasing gradually Wh or a, respectively. 
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The profiles of the metric functions of a typical EGB-AdS black hole solution corresponding to a = 0.5, 
r/i — 0.5 are presented on Figure 1 for a static (w^ — 0) and a rotating solution with — 1.8. One can see 
that the rotation leads to non constant values for hir^jr^ and &(r)//(r), and is particularly apparent on the 
function h(r) and its derivative. 

It is also natural to study how the profile of a rotating solution (e.g. with a given angular velocity w-h) 
is affected by the Gauss-Bonnet term. This is illustrated on Figure 2 where the profiles corresponding to 
a = 0.1 (dashed curves) -very close to the MPAdSs solution, and a — \ (solid lines) are superposed for 
r/i — 0.5, Wh = 0.5. One can see there that the r^— terms start dominating the profile of the metric functions 
/, 6, h very rapidly, which implies a small difference between different solutions for large enough r. However, 
the situation is different in the small— r regime (see also Figure 1). 

We also performed an analysis of the EGB-AdS solutions when varying the Gauss-Bonnet coupling 
constant a. In the limit a = 0, the MPAdSs black holes ()2.15|) are recoverecFl. The evolution of the 
parameters /'(r/j), b'{rh)^ h{rh) and J2, b2, characterizing the solutions is shown on Figure 3a as function 
of a. The corresponding physical quantities, as computed according to the relations in the previous Section, 
are reported on Figure 3b. 

We also varied the event horizon radius for a set of given a, Wh and found no evidence of a maximal 
value of rfi where the solutions could eventually terminate. The evolution of the solution data as a function 
of the event horizon radius is reported on Figure 4a, the mass-energy, angular momentum, entropy and 
Hawking temperature being plotted in Figure 4b. For small values of the numerical analysis is quite 
tedious and it strongly suggests that the derivatives of w{r) and h{r)/r'^ become infinite in the limit. 

Finally, although the numerics is more involved in this case, we constructed solutions with fixed a and rh 
but varying the horizon velocity wi^. Equivalently, this leads to a family of solution with varying the angular 
momentum J W4 since there is a one-one correspondence between and Wh- Similar to the a — case, 
for each set of solutions we observe two branches, extending up to a maximal value of Wf^ where they merge 
and end. The lower branch emerges from the static solution in the limit = 0. The maximal value of Wh 
depends on the horizon radius r^, the cosmological constant A, and the coupling constant a. 

5 Further remarks 

The main purpose of this paper was to present arguments for the existence of rotating black holes in d = 4 + 1 
EGB theory with negative cosmological constant. These configurations posses a regular horizon of spherical 
topology and have two equal-magnitude angular momenta, representing generalizations of a particular class 
of MPAdSs black holes. We also proposed to adapt the boundary counterterm formalism of [10] to d = 4 + 1 
EGB-A theory, computing in this way the mass-energy and angular momenta of solutions. The general 
relations in Section 2.1 apply also to other known solutions in EGB theory with negative cosmological 
constant and can easily be generalised for a positive sign of A. 

The solutions in this paper may provide a fertile ground for further study of rotating configurations in 
EGB theory. For example, their generalization to include the effects of an electromagnetic field is straight- 
forward. Also, in principle, by using the same techniques, there should be no difficulty to construct similar 
solutions in d = 2N + 1 dimensions with N > 2 equal magnitude angular momenta. An interesting problem 
here is to find the boundary counterterm expression in EGB— A theory for other values of d > 5. 

In the five dimensional case, one can also approach the general case of a black hole with two distinct 
angular momenta, by solving a set of partial differential equations with a dependence on (r, 6). The formalism 
proposed in Section 2 to compute the mass, angular momentum and entropy of AdS solutions should apply in 
the general case, too. Rotating topological black holes in EGB theory with an horizon of negative curvature 
are also likely to exist for A < 0. 

The study of the solutions discussed in this paper in an AdS/CFT context is an interesting open question. 
According to the AdS/CFT correspondence, the higher derivatives curvature terms can be viewed as the 
corrections of large N expansion of the boundary GET in the strong coupling limit. For the ansatz considered 
here, the boundary metric is not rotating and corresponds to a static Einstein universe in four dimensions. 

^'^With the particular values rjj = 0.5, ui^ = 0.5 that we have chosen to perform the numerical analysis, the parameters a, M 
of the MPAdSs solution correspond to a = 10/81 and M = 6561/48640. 
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Here it is interesting to note that, similar to the a = case, the stress-energy tensor for the dual theory 
defined in that background, as computed according to the standard prescription ^7\, is traceless. 

A detailed study of the d = A + 1 rotating black hole solutions in EGB theory together with a discussion 
of their asymptotically flat limit will be presented elsewhere. 
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Note added: When this work has been at final stages, the authors of '28^ appear to have succeeded in 
constructing a Kerr black hole in EGB theory in d = 4 + 1 dimensions. However, this result needs an 
independant confirmation. Afterwards, it would be interesting to extend the results in Section HI of |28| to 
the case of black holes with negative cosmological constant and two equal magnitude angular momenta. 
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